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1. Laplace- type differential operators 

The convential Laplacian A = div grad on the Euchdean space H" admits a natural 
generalization to a Riemannian manifold {M,g): 

A = g'^'VaVb. 

Here 

g = gabdx''dx\ (g-^) -.^ (g^,)-^ 

with respect to local coordinates , and V is the Levi-Civita derivative 

to g acting on scalar or tensor fields. Another far-reaching generalization is a Laplace-type 
linear differential operator of the form 

L = g'^^'DaDb + W (1) 

acting on the smooth sections of a vector bundle E over the manifold M. A covariant 
derivative D of these sections and a section W of the endomorphism bundle EndE to E 
enter the expression for L. 

An intrinsic definition, alternative to the coordinate-dependent definition (1), can be 
given [54]. Let the map adf for a smooth scalar field / e C°°(M) act on any linear 
operator L according to 

adfL := [L,f] = Lf-fL. 

(Here and in the following a scalar field is identified with the operator of multiplication 
by this scalar field.) A linear operator L acting on (p e C°°{E) is called a Laplace-type 
differential operator if 

{adffL^2\dff. 

It can be shown then that 

2 < (if, Dip >:= {adfL)ip - (A/)v? 
defines a covariant derivative D on C°° (E) and that 

Wif := Lcp - g'^^DaD^if 
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defines an endomorphism W e C°°{EndE). We use the notation 

<v,w >:= g°'^VaWij, |f p :=< v,v> 

for bundle- valued one- forms v = Vadx"', w = Wadx°'. 

Clearly, the class of Laplace-type operators L is form-invariant under diffeomorphisms 
when the objects g, D, W are carried along with This class is also form-invariant 
under gauge transformations 

L := A-^LA, 

where A G {Aut E) is a section of the automorphism bundle to E, and under conformal 
transformations 

L := e-^'^+^^fLe'^f, 

where / G C°°(M) and m := {d — 2)/2, d being the dimension of M. 
We have 

g = g^ D = A-^DA, W = A-^WA 
under a gauge transformation with A, and 

g = e'^fg, W = e-'^f {W + I ■ e'^^f Ae'^f) 

under a conformal transformation with /. Multiple covariant derivatives D of sections of 
E behave like multiple Levi-Civita derivatives of a scalar field under conformal transfor- 
mations. In particular, 

Dip = Dip for ip G C°^{E). 

For details cf. [54,33]. 

2. The heat kernel coefficients 

Some sequence of locally defined two-point quantities Hh — Hk{x, x') {k = 0, 1, 2, . . .) 
is associated to any Laplace- type operator L. 

Proposition 1. A formal series expression 

K = K{t;x,x') = (Ant)-"/^ exp f^tl^\ Y.Hk{x.x')t^ (2) 

^ ^ fc=0 

solves the heat equation 

dt 

with the initial condition 

K{+{);x,x') = 5{x,x') 

if and only if 

(P + //)iyo = 0, Hoix,x)=I, (3) 
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and 

(V + k + ii) Hk = LHk-u for A; = 1,2,.... (4) 

Here a{x, x') is the so called geodetic interval defined as one half the square of the distance 
along the geodesic between the points x and x', I> is a first-order linear differential operator 
given by 

T>(p =< da, Dip >, 

/ is the identity endomorphism, and = |(Acr — d). 

All the analysis in this paper is purely local. This means that we fix some small 
regular region of the manifold M and consider the points x and x' to lie inside this region. 

Proposition 2. There is a neighbourhood N of the diagonal of M x M and a unique smooth 
system of solutions Hk = Hk{x,x') in of the differential recursion system (4). These 
Hk — Hk{x,x') are called the Hadamard coefficients to L. 

Note that the quantities H]^ = Hk{x,x') are also widely known under the name 
Hadamard-Minakshisundaram-De Witt-Seeley (HMDS or HAMIDEW coefficients), ac- 
cording to the papers of these authors [36,42,43,21,22,51,52], or heat kernel coefficients, 
the last being used in the title of this paper. Hadamard [36] introduced them already in 
1923 for scalar operators L and established the essential properties. 

3. Methods for effective calculations 

We want to classify the methods for an effective calculation of the Taylor coeffi- 
cients [D'^Hk] {p-ik — 0,1,2,...) of the Hadamard coefficients or, which means a little 
more, of the diagonal values of their multiple covariant derivatives with respect to x, 
[Daj^Da^ • • • Da^H).] (n,k = 0,1,2,...). The square brackets [ ] applied to a two-point 
quantity means the restriction to the diagonal of M x M. Actually, we arrange the meth- 
ods to five groups. 

I. Taylor expansion in normal coordinates. Naturally, this is a method which 
comes first to the mind. It has been used with great skill by Giinther, Wiinsch, McLe- 
naghan and others [31-33,65,41]. 

II. Manifestly covariant methods. The Taylor expansion in normal coordinates 
can be reinterpreted in a coordinate-independent manner. So, if one works with covariant 
derivatives instead of partial derivatives it should be possible to replace the method I. by 
manifestly covariant formulas. In fact, Synge [60,61] and De Witt and Brehme [21,22] 
made the first steps in this direction. Their work can be continued in two ways to an 
algorithm for the Taylor coefficients [D'^Hk], which are given by the diagonal values of the 
symmetrized covariant derivatives of Hf,, or to an algorithm for the diagonal values of the 
{non- symmetrized) covariant derivatives. We will present such algorithms, developed by 
the authors, in the next sections. 

III. Invariant-theoretical methods. Let Rabcd denote the components of the 
Riemannian curvature of (M, g) and Fab be the components of the curvature of the bundle 
connection D, 

[Da,Db]<f = Fab<f fOT<feC°°{E). 
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The components of [D^Hk] or of the diagonal values of the multiple covariant derivatives 
of the iffc, [Da^Da^ ■ ■ -Da^Hk], exhibit a controlled behavior under diffeomorphisms and 
gauge transformations. This implies that they are just polynomials in Qab, g""^, Rabcd^ Fab^ 
"^aiRabcd, Da^Fab, Vai Vaj-Rabcd, Da^Da^Fah, ctc. Some careful analysis shows that the 
coefficients of these polynomials are determined by the functorial properties under several 
decompositions: Riemannian product {M, g)—{Mi^ gi) x {M2,g2), direct vector bundle 
sum E = El (B E2, tensor product of vector bundles E = Ei (E) E2, and splitting of the 
potential term W = Wi + W2- Gilkey [26-29] extensively applied the invariant-theoretical 
method as we call this analysis. Let us mention also the papers [65,33,25]. 

IV. Heat semigroup methods. If L is elliptic, i.e. the metric is positive definite, 
and if the manifold M is compact, then the linear operator exp(tL) is well defined and forms 
a semigroup, the so-called heat semigroup. (Let us mention that there is also a construction 
of the heat semigroup for non-compact M under certain additional conditions.) The linear 
operator exp(tL) has a kernel K = K{t; x, x'). The latter admits the formal series (2) as 
an asymptotic series as t ^ +0, cf. [42,43,18,28]. 

In quantum field theory there have been developed several perturbational techniques 
for handling a semigroup of the form exp(tL). Let us mention the Feynman path integrals, 
graph methods, the Dyson formula etc. Such methods can produce results for the kernel 
K = K{t;x,x') of exp(tL) and, hence, for the Hadamard coefficients. Osborn, Zuk, 
Nepomechie and others used this approach [44-48,69-72,23]. In recent papers of one of 
the authors (I.G.A.) the whole heat kernel diagonal K{t;x,x) (not only the Hk{x,x)) in 
low-energy approximation was constructed using the semigroup operator approach [8-14]. 

The semigroup exp(tL) has also a probabilistic interpretation: it describes some 
stochastic process. In particular, exp(tL) on scalar fields describes Brownian motion; 
the general exp(tL) belongs to some more complicated Markov process. The probabilistic 
approach gives further insights; let us only mention the Feynman-Kac formula and let us 
quote the papers [17,39,19,20]. 

V. PSEUDODIFFERENTIAL OPERATOR METHODS. Hadamard in his book [36] con- 
structed an asymptotic expansion for the Green function to L; he called it elementary 
solution. Let L be elliptic and replace it by L — XI, where the parameter A belongs to 
the resolvent set of L, i.e. A is not a spectral value of L. Hadamard's construction for 
L — XI can be reinterpreted in terms of pseudodifi^erential operators as a parametrix con- 
struction. The coefficients of the asymptotic series are just the Hadamard coefficients, up 
to numerical factors. Knowing this, one can produce results for the Hk {k = 0, 1, 2, . . .) 
by means of refined pseudo differential operator methods. Gilkey, Fulling and Kennedy, 
Gusynin, Obukhov and others did this succesfuUy [27,28,24,35,46,16,64]. 

Note that our systematics considers only general methods. Additionally to them, 
there are special methods for restricted classes of Riemannian manifolds, vector bundles 
and/or differential operators. Let us only mention the harmonic analysis on Lie groups 
or homogeneous manifolds: a spectral analysis of group-invariant differential operators by 
means of representation theory. 

4. An algorithm for the Taylor coefficients of the Hk 

The exposition in this section is due to the papers of one of the authors (I.G.A.) [2,3,6,7]. 
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Hq{x,x') = A^/'^{x,x')r{x,x'). 

where A(a;, x') is the Van-Vleck-Morette determinant and V{x, x') is the parallel displace- 
ment operator with respect to D along the geodesic. The modified coefficients 

Ofe •— — — -"0 -"fe 
fulfill the simplified differential recursion system [22] 

ao = /, {l + au = Mak-i (5) 

where 

M := -Hq^LHq. (6) 

We assume that there is a unique geodesic connecting the arguments x and x' of 
our two-point functions and consider all two-point quantities to be single-valued analytic 
functions. 

The formal solution of the recursion relations (5) has the form 

afc=^l + ip^ m(^1 + -^^V^ M---{1 + V)-^M 

Let us expand afc in a covariant Taylor series 

Ofe = ^ \n >< n\ak >, 

n>0 

where some compact notation is used. Namely, let grad'a denote the vector field with 
components 

(grad'a)'^ := grad'a <^ ■ ■ • <S> grad'a 
be the n-fold (symmetric) tensor product and 

(-1)'' 

\n >:= ^— ^{grad'a)''. 

Let further the linear functionals < n| (dual to the system of functions \n >) be 
defined by 

< n\(fi >:= [D'^(f] 
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for a two-point field (p, where D"^ = DD ■ ■ ■ D is the multiple symmetric covariant derivative 
with respect to x, and the square brackets [ ] denote, as before, the restriction to the 
diagonal. Then |n >< n\(fi > denotes the inner product 

\n >< n\cp >= {\n >)'^i-^^"(< n\cp >)^i...^„. 

From the remarkable eigenvalue formula 

Pjn >= n\n > 

one can obtain an expression for the Taylor coefficients of the a^, namely 

ni,-",nfc_i>0 

X < n\M\nk-i >< nk-i\M\nk-2 > ■ ■ ■ < ni\M\0 >, 
where < m\M\n > are the 'matrix elements' of the operator M (6) 



< mlMln >- 



D^'M^—^igrad'ay 



(7) 



The matrix elements < m\M\n > are tensors with components M'^i - ''", „ which are 
symmetric both in upper and lower indices. The inner product of the matrix elements is 
again such a tensor with components 

(<n\M\k>< k\M\m>r'-''- ^ = M'''-''\ . M^'-\ „ . 

These tensors are expressible in terms of diagonal values of derivatives of basic two-point 
quantities. The components of these matrix elements were found in [2,3,6,7]. 

In order to present explicit formulas for these matrix elements we use a calculus 
of matrix- valued, vector- valued and EndE-valued symmetric differential forms. So we 
introduce first matrix- valued symmetric forms -f^(n) = (-^^^^(n))" 

K\^n) = V(^, • • • V^„_,i?-^_^|^|^„)dx^^ V ■ ■ • V dx^-, 

where V denotes the symmetric tensor product of symmetric forms. Then we define the 
matrix- valued symmetric forms 7(„) = (7"^(n))' ^(n) = (^"/3(n)) ^nd {X^"^^^) by 

l<A;<[f] ni,---,»Xfc>2 
'T'l H l~'^fc — 

where 

iV(ni , • • • , nfe) = (m - 2) ! • • • (n^ - 2)! 

(n- 1)! 

X ni(ni + l)(ni + n2)(ni + n2 + 1) • • • (ni H h nfe_i)(ni H h Uk-i + 1), 
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V{n) = - 2^ 2^ , 7(n.) V---V7(^,) V7(„,), 

l<fe<[f] "l.---."fc>2 

'T-l H 



{n-ky 



0<k<n 

and scalar symmetric forms (^n) 

),! 



Til H \-ni^=7i 



Next, let us introduce the vector-valued (and E'nd-E-valued) symmetric forms = 
(-^M(n)): 

and = (^^(n)) 



l<fe<n-2 ^ ^ 



and £^n(i£^-valued symmetric forms VF(n) 

Using the introduced quantities we define finally some more symmetric forms 

X-l^n)=X''l,,...,J^'''^---^dx>^-, 

and 

Aaf3{n) = Aa,pt,,...,,Jx>'' V • • • V dx>'- , Ca(3(n) = Ca^Mi-M.^^a;'^' V • • • V dx"" . 

where X'^^^ ^i...p„ and Cmi-.-a^™ are the components of the forms X^^^^y and 

C{n)- The components of the matrix elements < m\M\n > are given then by 

<m\M\n>=0 for n>m + 2 and n = m+l, 

n — 1) (M1---/^k-2 Mn-l-"Mm)' 
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where X'^^^...^^, ^^^^^...^^^ and Z^^...^^ are the components of the forms X^^ {n) 

y\n) = -I X^\in) + 2 E (^) X^-^,^ V , 

0<fc<n ^ ^ 

+ X] ^"(fc) ^ (-^a^(n-fc) + I CaPin-k)) + ^^''a(fc) ^ (-^/3(m) + C^(m)) [ 

0<fe<n \ ^ 

+ X! /cimlfn- fc-m)!^"?^) ^ ^ ^^(n-fc-m) - Ca(m) V C/3(n-/c-m)) • 

5. An algorithm for the derivatives of the Hk. 

The following is due to the papers of the other author (R.Sch.) [53-58]. Let us begin 
with a new tensor notation: indices at a tensor symbol shall not denote the components 
but the valences, that means the entries of the tensor taken as a multilinear functional. 
Positive integers 1,2,3... are preferred indices for the valences. Thus, u\2...p denotes a 
covariant p-tensor, v}'^" '^ denotes a contravariant p-tensor, Wi2...p + Vi2...p is the sum of 
two tensors, 'Ui2...p^^p+i,...,p+q is a tensor product, the natural action of a permutation tt of 
1, 2, . . . , p on a p— tensor is 

T^Ul2...p = 'U7r(l)7r(2)...7r(p)- 

More generally, if 11 is a subset of the symmetric group Sp (of permutations of 1, 2, ... 
then 

nWl2...p ^ TTU12...P- 

Tren 

We are interested here in special permutations: a qf-shuffle of 1, 2, . . . ,p is a tt e 5'p 
such that 

7r(l) < 7r(2) < ■ ■ • < 7r(g) and 7r{q + 1) < 7r{q + 2) < ■ ■ ■ < n{p). 
The action of the set of all q'-shuffles on a covariant tensor tti2...p is denoted by 

U l2...q q+l...p - 

For example: 

Ul23 = U123 + U213 + ^312, 
Ul2 34 = 1*1234 + ■J^1324 + ■J^1423 + ■J^2314 + ■J^2413 + ■J^3412- 

Let us define tensors Sai2...p ip > 3) and i?n(i i?- valued tensors Mi2...p {p > 2) by 

'S'al23 := Ral23, 'S'al234 := -Ral23;4 
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p-3 

Sal2...p '■— -Rol23; 45---p + Ralb2; 3---q S^q+l...p, 

-^12 Fi2, Mi23 :— -?^12;3) 

P-3 

M12...P := i^l2; 34...p + ^ i^al; 23...g 'S'°g+i...p, 

where indices after a semicolon express covariant derivatives. We define also 

Sal2...p '■= [C;al2...p] + Sal2...p 

and 

1^12. ..p := [/^;12...p] - ^12. ..p- 

The diagonal values of the derivatives of the Hadamard coefficients are then recursively 
given by 

{p+k)[Hk;12...p] = [{LHk-l);12...p)] 

p p 

- ^ "^ 12...p [-^fc; q+l,...p ] - ^ s \2...q [Hk;a q+l...p ], 
q=2 q=3 

where 

{LHk);12...p = ^^fc;"al2...p + ( ) ^; 12...g -^fc; g+l...p - 

6. Discussion 

Let us add some historical remarks on the explicit calculation of the diagonal values of 
the Hadamard coefficients and their Taylor coefficients. Already in the thirties Heisenberg 
and Euler [37] and Mathisson [40] have calculated [Hi], [DHi], [D^Hi] for Laplace-type 
operators in flat space, L = g'^^DaDf, + g"'^ = const. For scalar operators in curved 
space let us mention Holder [38] for [Hi], Giinther [31] for [Hi], [H2], Sakai [50] for [H^], 
Amsterdamski, Berkin and O'Connor [1] for [-^4], Giinther [32] for [DHi] and [D^Hi] and 
Wiinsch [66-68] for [DH2] and [D'^H2]. For non-scalar operators in flat space we mention 
[62] for [H4] and [63] for [Hr^]. 

Results for the general Laplacc-typc operator have been obtained by Dc Witt 
[22]: [Hi], [H2], Gilkey [26]: [Hi],[H2],[Hs], one of the authors (I.G.A.) [2,3,6,7]: 
[Hi],[H2],[H3],[H4], the other author (R.Sch.) [53,54]: [DHi],[D^Hi],[D^Hi], Wiinsch 
[65]: [D'^Hi], and others. 

The Hadamard coefficients H^ or rather the one-point quantities [D^Hj^] or 
[Da^Da2 ■ ■ ■ Da^Hk] derived from them have many appearances or applications: 

— Huygens' principle for hyperbolic L [40,31-33,41,53-55,65-68], 

— Spectral geometry for elliptic L and compact Riemannian (M, g) with positive definite 
metric g [42,43,18], 
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— "Heat kernel proofs" of the index theorems [34,28], 

— Zeta-function regularization, 

— Regularized energy-momentum tensor, 

— Korteveg-De Vries hierarchy [56,15]. 

It should be noted that there are some other problems in mathematical physics that 
require a similar technique for investigation of transport equations. These are: 

— expansion of the metric and other quantities in normal coordinates [33] , 

— harmonic manifolds [49] and harmonic differential operators [58], 

— volume problems in the sense of [30] (to read geometric information from the volume 
of geodesic balls, truncated light cones, ...) [30,57], 

— Brownian motion. 
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